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Abstract 

By using the random matrix approach and generahzed Blue's function representation 
we solve analytically the model of an effective two-level system coupled to a noise 
reservoir. We show that calculated spectral properties of the system are in agreement 
with the numerically simulated results. We outline possible applications of the model 
in the field of condensed phase reactions. 



1 Introduction 



The theory of open quantum systems has been found to play a key role in many ar- 
eas of theoretical physics |l|]. The growing interest focuses on quantum dynamics in 
condensed phase systems, to which the most widely used approach is the reduced- 
density-matrix theory [H, |, ^, H. While a microscopic treatment of dissipative motion 
involves hermitean Hamiltonians, use of non-hermitean Hamiltonians serves as a con- 
venient way to give a reduced description where the hermitean part of the Hamiltonian 
refers to free undamped dynamics and the remainder describes a damping imposed on 
the system by some external noise source ("heat bath"). 

Statistical properties of complex systems can be successfully investigated within the 
framework of the random-matrix theory which turned out to be quite general and 
a powerful phenomenological approach to a description of various phenomena such as 
quantum chaos [0, complex nuclei 0], chaotic scattering and mesoscopic physics [Q. 
Aspects of vastly different physical situations such as electron localization phenomena 
in disordered conductors and semiconductors B, disordered quantum wires [|T0| and 



quantum Hall effect can be described in the language of the random matrix theory. 
In all the realms mentioned above, the Hamiltonian of the system is rather intricate to 
be handled or simply unknown. In such cases the integration of the exact equations is 
replaced by the study of the joint distribution function of the matrix elements of the 
Hamiltonian -P(H). If there is no preferential basis in the space of matrix elements, 
so that the system is as random as possible and equal weights are given to all kinds of 
interactions, one has to require the probability P(H)(iH to be invariant under similarity 
transformations H — > R~^liR with R being orthogonal, unitary or a symplectic N x N 
matrix reflecting the fundamental symmetry of the underlying Hamiltonian. 

The natural way of addressing the problems of randomness coupled to various sources 



is to use technique of the "free random variables" [T^, |13|, 0. This method pro- 
vides an elegant way of "linearizing" the process of determining the average eigenvalue 
distributions for convolutions representing an analogue of the logarithm of Fourier 
transformation of the usual convolutions. Recently, the generalization of the "addition 
law" for hermitean random matrices to the non-hermitean case has been derived [15|. 

In section 2, we introduce, after Zee [|l^] and Janik et al. ||15| basic concepts of the 
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"addition laws" , for both hermitean and non-hermitean ensembles. 

In section 3, we apply the addition law to the dissipative system of a two- level deter- 
ministic Hamiltonian coupled to a random noise. The effective Hamiltonian is obtained 
by the standard Wigner-Weisskopf reduction [0, |1^, |18[ by partition. In this method 
the total Hilbert space is divided into two subspaces and the Hamiltonian is integrated 
over the elements of one of them, eventually mapping the eigenvalue problem of high 
dimension onto a lower one. Models of such type are widely used in e.g. condensed 



phase dynamics [0, ^ where one is typically interested in the detailed dynamics of 
only a small part of the overall system, so one partitions the total system into the 
subsystem of interest and a bath: 

H = Hgys + Hbath + V (1.1) 

Here Hgyg and Hhath are the Hamiltonians of the isolated subsystem and the bath; 
V describes their interaction. The appropriate scheme for partitioning is not always 
obvious [0; successful choice minimizes the coupling V so that its effects can be studied 
perturbatively. After the separation, the state of the system is described by its reduced 
density matrix which is obtained from the complete system-plus-bath density by taking 
the trace over the bath variables. When the relaxation of the bath is rapid compared 
to the evolution of the system, the bath can be assumed to be equilibrated at all times. 
Treating the interaction of the system and bath perturbatively to the second order, 
the evolution of the reduced density matrix can be then recasted in the form of the 
dissipative Liouville-von-Neumann equation 0|. 

Similar ideology is used e.g. to study the electron or wave passage in disordered 
medium. In the field of chemical physics the condensed phase reactions, such as electron 
and proton transfers are of primary interest and are commonly investigated in terms 



of a dissipative two- level system, often referred to as spin-boson problem [19, pi[. 



We solve the analogue of this problem by investigating analytically an effective two-level 
system coupled to a Gaussian dissipative term. We demonstrate the phase transition 
in such a system corresponding to the structural change in the average spectral dis- 
tribution of the eigenvalues. We present also the results of the numerical simulations, 
confirming our analytical results. Finally, section 5 outlines the possible applications 
and generalizations of the considered model. 



2 



2 Addition laws for random ensembles. 



The generic problem of finding the Green's function of a sum of two independent 
hermitean matrices (here Hamiltonians) of dimension N 

01 + 02 (2.2) 



has been considered in the seminal work by Voiculescu |]T2|, and recently popularized 
by Zee 0. 

The main idea of Voiculescu was to find an elegant way of "linearizing" the process of 
determining the average eigenvalue distribution for the convolution of non-commuting 
operators, hence finding an analogue of the logarithm of Fourier transformation for the 
usual convolutions. As an example, let us consider two hermitean operators in matrix 
representation, (|2.2|) , populated from the statistical ensemble P(0i,02), provided that 
in the large N limit P(0i,02) factorizes into -Pi(0i) • -P2(02)- Then the resolvent 

G{z) := / c?0irf02P(0i, 02)Tr (2.3) 

N J 2-01-02 

could be expressed (at least a priori) knowing only the individual resolvents 

G,(2) = 4 /#^^*(0.)Tr^- i = l,2. (2.4) 

i\ J z — (pi 

The physical interpretation of the addition law comes from the observation that 

G(z) = — Tr ( = (2.5) 

^ ^ N \z-(t)l z-T.{z) ^ ' 

where 

(...)=/ P(0) . . . [rf0] (2.6) 



and S(2;) is nothing but the self-energy. Through his diagrammatic analysis. Zee |T3 
introduced the "Blue's function" that is just the functional inverse of the resolvent 

B{G{z)) = z (2.7) 

and satisfies the additivity law 

B,+2iz) = B,{z) + B2iz)-- (2.8) 
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The operational algorithm is extremely simple when using the Blue's function approach: 
First, having both ensembles, construct the two resolvents Gi and 6*2- Second, invert 
them functionally, resulting in i?2- Third, add them using the "addition law" ( ^.S] ). 
Fourth, invert functionally the -B1+2, getting the Green function for the sum G1+2 = G. 
Fifth, get the average spectral distribution from G, using the standard construction 

g{e) = -- lim Im G{z)U=,+a- (2.9) 
vr A^o 

The generalization for the nonhermitean case amounts to considering the isospin 
Green functions which now form a 2 by 2 matrix in isospin space, indexed by z and z. 
The generalized Blue's function ||15| is now a matrix valued function of a 2 x 2 matrix 
variable defined by 

m = ^ = ( I) ^2.10) 

The Green's function is found from the matrix version of (12T 



(0 lj=B^ig) + B,ig)-^ (2.11) 

In general, two kinds of the solutions to this matrix equations are possible. The first, 
when Q is diagonal, reproduces the hermitean construction modulo the anti-analytical 
copy in z, totally decoupled from the 2;-solution. Therefore this one reproduces the 
hermitean (or antihermitean case) only. The relevant solution is the one when the 
off-diagonal element of the Green function is non-zero. In this case, the resolvent is 
nonanalytic, being the function of both variables z,z. The support of the eigenvalues 
forms the two-dimensional "islands" of nonanalyticity, (opposite to the hermitean case, 
when only cuts constitute the support of real eigenvalues) and the density distribution 
could be read-out from the 2-dimensional "Gauss law" 



divE = 47r^(2,^) (2.12) 

In this two-dimensional electrostatics divergence means differentiation with respect to 
z, and electric field has two components, = 2Re^2^, Ey = —2lmQzz- 

Note that the non-analytical solution is a consequence of the "spontaneous breakdown 
of analyticity" |]I5|. Naively, the resolvent ( p.5| ) is the function of z variable only. 



But the true ground state is not respecting holomorphic separability into z and z 
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copies, resulting in the spontaneous breakdown of this "symmetry", by forming the 
mixed "condensates" Qzz and therefore breaking the analyticity (non-analyticity) of 
the solution[|. Below we demonstrate this scenario in the case of Wigner-Weisskopf 
reduction. 



3 Dissipative two-level system. 

We consider the Hamiltonian with the Hilbert space spanned by discrete states \k) 
and M continua \xn{E)): 

N M „ 

n = Y.Ekm{i^k\ + E / dEE\Xn){Xn\ 

k=l n=l-' 

M N 

+7^/' E E / dE{V-{E)\^,){xn\ + h.c.) . (3.13) 

n=l k=l •' 

The bound continuum coupling is characterized by the coupling constant 7 and the 
energy dependent matrix V{}{E). Note that in this model there is no direct coupling (in- 
teractions) between continuum channels n,n' . After eliminating the continua \xn{E)), 
the effective Hamiltonian H^ff takes the form 

H,sf = n-iiVV\ (3.14) 

where Ti. is N x N deterministic Hamiltonian with the bound-state energies ei, . . . ,6^ 
and the dissipative part is obtained via so-called Mar kovianQ limit. The operators VJ? = 
{ipk\V\Xn) stand for the coupling of the bound-state eigenvectors ipk {k = 1, . . . , N) cor- 
responding to eigenvalues ei, . . . ,eN with the continuum states Xn {n = 1, . . . , M) and 
account for the depletion of the bound state subspace due to the leakage of the proba- 
bility mass into the continua. The eigenvalues e = x — iy of the effective Hamiltonian 
Heff are in general complex, their imaginary parts must be positive for damped systems 
and give the lifetimes of the corresponding eigenvalues proportional to 1/y. 

*In other words, the free energy associated with the non-analytical solution lies below the one 
related to the analytical resolvent in the domain of the complex plane. 

^The Markovian property assumes exponential relaxation of the response of the continuum when 
it acts on the test system and rephrased in thermodynamic terms, would correspond to the high tem- 
perature Debye relaxation bath. In the above model, Markovian approximation should be understood 
as the condition for almost instantaneous decay of the kernel in the evolution (Schrodinger) equation 
satisfied by the bound-space state vectors Upk)- The condition is met if the energy dependence of the 
matrix elements (E) is negligible [|[ |[ || . 
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The model Hamiltonian of that type ( |3.14| ) is widely used in quantum chaotic scat- 
tering problems ^ but has been also introduced to study unimolecular dissociation 
of selectively excited polyatomic molecules E^, HM. For most dissociating molecules 
above the tunneling regime (when the excitation energy is well above the dissociation 
threshold) the coupling between the molecular complex and the continuum is strong 
and the density of strongly mixed molecular states is large, so that the unimolecular 
decay is characterized by overlapping resonance states. 

We will solve the system in the regime when both M and are large, but m = M/N 
is fixed (which implies strongly overlapping resonances). To simulate the effect of the 
"noise" represented by statistical properties of the operators V, V"^ we will assume that 
the N by M matrices V are random. For simplicity we will choose randomness as 
Gaussian 0, thus working within the Gaussian unitary ensemble (GUE). Furthermore, 
we will investigate a two-level system, with deterministic eigenvalues ei = . . . = eN/2 = 
e and eN/2 = ■ ■ ■ = = -e. 

To apply the addition law ( |2.11| ) , we need only the functional form of the corresponding 
Blue's functions. For the deterministic part, the Green function is by definition 



2\Z-e Z 

where z has been replaced by the matrix Z = diag(z, z)). Therefore the relevant Blue's 
function is given by the equation 

Z = i {{Bn{Z) - e)-' + (BniZ) + e)-') . (3.16) 

The Blue's function for the noise term i'-yVV'^ can be determined by use of diagrammatic 
techniques and has been first derived in [|T^: 

Bvvt =mil-TZ)-'T + ^. (3.17) 

Z 

where F is the coupling matrix 

r=(~r I). (3.18) 



*By use of the central limit theorem, Gaussianity is a legitimate choice for the infinite normalized 
sum of independent fluctuating contributions. It can be shown also that the distribution for normalized 
sum of matrices 0sum = '^^i'Pi' each sampled from a distribution Pi{4>i) follows (l3[ Q the Gauss 
theorem of probability. Using measures different from the Gaussian one do not change the addition 
law discussed in this section. 
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Note that all above equations have been written in the matrix form, in agreement with 
the discussion of the previous section. 



We can further use the generalized addition formula for the total system, result- 



ing m 

B{Z) = Bd{Z) + m(l - TZy^T . 
Inverting that relation {i.e. replacing Z hy Q), we get 

Z = BD{g) + m{l-Tg)-'T . 
By similar replacement in equation (|3.16|) we get 

G = \[{Bo{G)-e)-' + {Bn{G) + er^ 



(3.19) 



(3.20) 



(3.21) 



Calculating Bd{G) from ( ^.2(JD and plugging it into ( ^.211 ) we arrive at the final equation 
for the Green function 



g = U{z- m(i - rG)-'r - e)-' + {Z- m(l - tgVt + t] 



-1 



(3.22) 



The matrix equation ( |3.22| ) is a central formula to the paper from which all further 
discussion of spectral properties of the system will be derived. We solve this equation 
using the explicit parameterization of the matrix G with the unknown entries 



a b 
d c 



(3.23) 



As already mentioned in the previous section, the solutions belong to two classes: 



• Holomorphic (analytical) case. 

This yields a trivial solution to the problem. In this case, we seek solutions with an 
imposed condition 6 = (i = 0. The above matrix equation splits then into two mirror 
copies of algebraic equations for z and respectively. The equation for a = Gzz is of 
Cardano type 

a%e^6^ - + a\6^z + 26z^ - 2m6^z - 2e^6) 
+a(m5^ + 2'm5z — m^5^ — 25z — z^ + e^) + z — m5 = 

5 = -i-i (3.24) 
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with the Blue's function 

Biz) 



m6 1 

+ 7^ + 



l-6z 2z 



1 

4^2 



1/2 



(3.25) 



Solution to those equations constitute nothing but a complement to the dashed domains 
of density g{z,z) pictured in Fig. 1. We reject it as a non-physical one. The proper 
solution to our problem follows a non-analytical case which we are going to discuss 
now. 




1 
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r ^^^^^ ' ^ 1 
1 1 1 



Figure 1: The support of the eigenvalues of the two-state system for different coupling 
at fixed filling ratio m = 0.25. The dots correspond to the numerical simulation of 
the Hamiltonian (|3.14|) , while the solid lines come from the analytical result ( p. 29 ). 
Numerical simulations have been performed for 100 matrices of size = 100 at 7 = 
1,7 = 3 (the upper row of drawings), and for 30 matrices with = 400 at 7 = 6 
and 7 = 10 (the lower row). The density of the dots is related to the the spectral 
density (^1281) . 



• Non-holomorphic (non-analytical) case. 

This case corresponds to a situation when, despite the resolvent naively is the function 
of the z variable only, the support of the spectrum depends on both z and z, forming 
therefore two-dimensional islands (in contrast to one-dimensional cuts of the analytical 
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solution) in the complex plane. The solution of the matrix equation now is given by 

"/I Ix^ — m — r 



X 



X 



a 

c 



+ 



+ 



7?/ 2(x^ ~ e^) 2 V7 7 y 



+ 



and b = d, where b satisfies the equation 

9 Tn I 1 9 1 ^ / \ 
b= — + |ar + ^ + -(a - a) . 

7Z/ 7 7 

The spectral density obtained from Gauss law ( p.l2| ) reads 



^ An 



All 



x"^ — 



2/^ 



+ 



m — 1 



2/^ 



(3.26) 



(3.27) 



(3.28) 



The shape of the envelopes (i.e. boundaries of 2-dimensional regions of spectral den- 
sities in the complex plane) is given by the condition b = Qzz = which is equivalent 
to 



X 



X 

T 



2 

ly 



1 \2 I 


[i(lH 


x^ -e^\ 


m — 1" 


2 


x^-e^J ~^ 4 


L7 V 


y' ) 


+ 

y J 





(3.29) 



The spectrum undergoes two various phase transitions (cf. Fig. p. For very low cou- 
pling, the energy islands do not connect the two states and only after reaching the 
critical value of the coupling 7, a single island of the eigenvalues is being created (see 
the upper row of drawings in Fig. |^). We note, that the global energy scale e may be 
eliminated in (|3.29| ) by introducing the dimensionless variables x' = x/e, y' = y/e and 
7' = 7/e- 

The second critical point (see the lower row of drawings in Fig. |I]) corresponds to the 
situation when one island of non-analyticity splits into two, a low lying broad island 
and a narrow island of short-living states. The critical values of the coupling constant 
can be determined from the study of the discriminant of the boundary equation ( |3.29| ) 
at a; = 0. The zero value of discriminant appears if distinct roots of the equation 
coincide, thus it is suggestive of occurrence of phase transitions discussed in the text. 
The condition is described by a second order equation in 7^, 



1 - mf-ft + (m^ - 20m - 8)e^-f'i + 16e 



0. 



(3.30) 



The solution 7* scales with the energy scale e, and the critical dimensionless coupling 
7*/e is universal. Figure ^ shows the two critical couplings as the function of the 
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filling ratio m. The lower value 7c, corresponding to the connection of the two energy 
states, starts at value 2 and goes to 4/3^3 as m 1. The other branch of critical 
coupling parameter 7^, describing the splitting of the broad low lying island and the 
narrow, short-living region, starts at the same value as 7c and goes to infinity as 
m — i> 1. For small values of m (m ^ 0), appearance of bridging states between the two 
"deterministic" states results in the immediate splitting of the cloud of the resonance 
states. Distribution of points within the clouds follows the ratio (1 — m) : m between 
the low lying broad cloud and the short-living upper cloud, respectively. 




' 1 ' ' 1 1 1 I ' 1 1 ' 1 

0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

m m 



Figure 2: The dimensionless critical couplings of the two-state model as functions of 
the filling ratio m. The left figure presents the critical coupling 7c at which the corridor 
between the two state is built up, the right figure displays the critical coupling 7^ at 
which the cloud of resonance states splits. 



By integrating g{x, y) over the support of x we obtain the distribution of resonance 
widths in the model. 

Using the symmetry of the spectra we arrive at 

0iy) = Jdx gix,y) = R{xu,y) - Rixi,y) , (3.31) 
support 

, 1 ( 2x 2x f x^ (1 — m)x X 1 

R{^,y) = ^\ — + — [^-^]-- — ^ + -2 — 2)' 

zTT [ 7?/ 72/'^ y 3 J e"^ — a;^ J 

where x^ and xi are the two positive end points of the support (for a continuous 
support = 0). The resulting density can be used to determine average decay rate in 
the system, and is shown in Fig. |^. The bound state component \ip) evolves in time 
according to 

|^(t)) = e-'^^^/*/'^|^(0)) (3.32) 
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Figure 3: The integrated spectral density, g{y) describing the density of the resonance 
widths in the system governed by -f^e//- Sohd hne stands for the analytical prediction, 
the dashed line represents numerical estimate obtained by use of an ensembles of 30 
random matrices with size = 400. 



and its effective decay rate can be defined as the logarithm derivative of its norm 



(3.33) 



If the initial state has nonzero overlap with only one eigenstate of H(,ff, the above 
definition simplifies to /cg// = 2?/j//i. In general case when resonances overlap, the 
average survival probability over the ensemble is given by 



Pit) 



Q{y) exp{-yt)dy 



(3.34) 



and is related to the distribution of kinetic rates g{y) by the Laplace transform. The 
result is displayed in Fig. ^ showing the decay with some average decay rate for low 
couplings and the appearance of roughly two distinct decay rates for large couplings. 
The time scale is of the order of a femtosecond for e = leV, and scales linearly with 
energy. (Note, that within the paper the mean-square dispersion of elements of V 
matrix has been put equal 1.) 



By integrating the spectral density (|3.28|) along the width coordinates one can obtain 
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Figure 4: The decay of the survival probabihty ( p.34 ) for various values of the coupling 



constant 7. Nonexponential decay becomes transparent for increasing values of 7. The 
solid line is the Laplace transform of the analytical width distribution ( |3.32|) . The dots 
correspond to the Laplace transform of the density obtained numerically with a set of 
30 random matrices of size = 400. 

the total width assigned to a given energy. The resulting density reads 

six) = Jdy g{x,y) = J2 R^i^^V'ii+i) - Rx{x,y2i+2) , (3.35) 
support *=i'2 

1 { — 1 — m + ln|?/|l 

R{^,y) = t{ — + ^71 — ^ + 

47r y 7?/^ y (e^ — x'^Y g J 

with (1/1-1/2) and {yz^yAj the supports of the integration. The inverse of that quantity 
gives information about the decay time of the state, labeled with the given energy. 
The distribution of the total width as the function of the energy at different couplings 
is presented in Fig. ^. The situation depicted in Fig. ^ strongly resembles "noise 
induced transition" when coupling to a fluctuating source induces new order in the 
system. For low values of 7, distribution q{x) has maxima around eigenvalues of the 
deterministic Hamiltonian; that bimodality of q{x) changes with increasing value of 7 
which introduces new maximum around x = 0, strongly populated for high values of 
the coupling constant. 
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Figure 5: The integrated spectral density, g{x) describing the total width of a given 
energy state. Solid line stands for the analytical prediction, the dashed line represents 
numerical estimate obtained by use of an ensembles of 30 random matrices of size 
N = 400. 

4 Collectivization of the spectra 



At large coupling constant the spectrum shows the structural change. Two groups of 
distinct eigenvalues appear. The first, saturating the most of the total width corre- 
sponds to the direct resonances. The second group corresponds to long-living states 
and appears only, if the coupling is sufficiently strong. Such a "reorganization" of the 
spectrum is usually observed [^, ^ in systems with N degenerate bound states 

coupled to M continua. The system develops then M decaying states and N — M sta- 
ble states whose population becomes trapped. The similar phenomenon of formation 
of a short-lived coherently decaying state in a system with almost degenerate levels 
is known in quantum optics (Dicke superradiance |29|), nuclear physics (giant dipole 
resonances, collective isobar- hole states in scattering at intermediate energies 0]) and 



dissipative spin systems [^. An interesting aspect of the phenomenon brings also a 
semiclassical description of collectivization discussed in papers of Gaspard and Rice 



3^ . Multiple hoppings lead to very long trajectories (the delay time is large) satu- 



rating the Gutzwiller formula and causing breakdown of the ergodicity. The physical 
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mechanism responsible for the appearance of the above mentioned structural change in 
the spectra is caused by the presence of zero modes in the matrix VV'^ . If the coupling 
is strong, the zero modes are becoming delocalized, corresponding to the appearance of 
the coherent, collective state. In other words, for increasing coupling 7, the probability 
of hopping from one zero mode to another (or effective "overlap" between the zero- 
mode states) is increasing, and the long strings due to multiple hopping can appear. 
This mechanism is somewhat analogous to the Mott's conductor-insulator phase tran- 
sition, when the conductance appears due to non-zero probability of the hopping of the 
electron between the various Fermi levels, so that the electrons are no longer localized. 
Another physical example corresponds to the picture of spontaneous breakdown of the 
chiral symmetry in quantum chromodynamics [^], when chiral condensate appears as 
a coherent effect of the delocalization of the zero modes of gluonic excitations. 



5 Conclusions 

Our analysis of the two-level system coupled to a random noise shows the advantage 



of using the Blue's function approach [|T3 to study spectral properties of dynamical 



quantum systems. By use of the generalized Blue's function |T5[, we are able to con- 
struct easily the Green function relevant for the system and from the last one we can 
read off spectral properties of the effective Hamiltonian. 

Characteristic features of the distribution function of the eigenvalues related to the 
effective Hamiltonian depend strongly on the intensity of the coupling between the 
deterministic and random parts of the Hamiltonian. At low values of coupling constant, 
we observe some dispersion of eigenvalues density around deterministic energies of the 
two-level system. Beyond the first critical value of 7 = 7c, a common region connecting 
"dispersed deterministic eigenvalues" is created and builds up with further increasing 
7. At very high values of coupling, dissipativity entering the system through the 
term VV"^ leads to reorganization of eigenvalues. We observe "collectivization" of 
widths to a bridging, narrow region which is populated by low energy states. Its 
presence suggests then formation of long-living states, well separated from the cloud 
of resonances responsible for a rapid decay. Despite its simplicity, the model presented 
here can cover vast number of various applications. We believe that one of direct 
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use of the analysis is modeling unimolecular decay processes in polyatomic molecules 
and electron tunneling processes in biological media. Studies of an effective tunneling 
matrix in biomolecules require knowledge of the structure of off-diagonal elements of 
the generalized Green function discussed above. Analysis of the kinetic rate in such 
models is based on the form of two-particle (two-time) correlation function which, for 
the model Hamiltonian of " deterministic plus random" dynamics, will be discussed in 
the forthcoming paper. 
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